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ABSTRACT 


The  existence  of  the  two  classes  of  thermal  neutrons  about  the 
Bragg  cut-off  in  a  crystalline  structure  has  been  experimentally 
verified  by  the  use  of  dysprosium-aluminum  wire  detector.   The  space 

dependent  neutron  flux  in  a  24x24x48  inch  rectangular  parallelopiped 

3 

graphite  pile,  GLC  grade  HC  and  density  1.60  gm/cm  ,  has  been  investi- 
gated at  295CK  and  481°K.   At  each  temperature  a  distinct  "knuckle"  or 
change  of  neutron  distribution  is  readily  apparent,  although  this 
degree  of  change  decreases  with  increasing  temperature  because  the 
neutron  cross  section  below  the  Bragg  cut-off  increases  with  increasing 
temperature. 
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I.   INTRODUCTION 

A  distinct  characteristic  of  the  neutron  scattering  cross  section 
is  known  to  exist  in  crystalline  materials  such  as  beryllium  and 
carbon.   The  cross  section  above  an  energy  level  corresponding  to  the 
Bragg  cut-off  is  relatively  constant.   However,  below  this  energy  level 
the  cross  section  experiences  a  sharp  drop,  dividing  the  thermal 
neutrons  into  what  may  be  called  the  cold  neutron  and  thermal  neutron 
groups.   Additionally,  within  the  cold  neutron  group,  the  cross  section 
varies  with  the  temperature  of  the  material  for  any  specific  neutron 
energy.   The  existence  of  these  groupings  is  graphically  illustrated  by 
the  variation  of  diffusion  length  across  the  thermal  neutron  energy 
spectrum.   The  temperature  dependence  of  the  cold  neutron  group  has 
been  studied  analytically  and  calculations  made  to  determine  the 
diffusion  lengths  in  graphite  at  temperatures  of  300°,  200°,  150°,  and 
100°K  by  Ahmed  [2]. 

It  was  felt  that  experimental  verification  of  this  theoretical 
temperature  dependence  was  in  order.  A  thermal  pile  was  constructed 
with  provisions  for  raising  the  graphite  temperature  up  to  650°K. 
Although  Ahmed's  calculations  were  based  on  below  room  temperature 
effects,  the  same  dependency  exists  for  the  higher  temperatures,  as 
can  be  seen  in  Table  1  which  includes  values  calculated  by  Trikha, 
et  al.  [21], 


TABLE  1 

TEMPERATURE 

(K°) 

THERMAL  GROUP  L 

COLD  GROUP  L 

REFERENCE 

100 

225.21 

141.75 

21 

150 

92.10 

45.00 

2 

200 

62.80 

43.10 

2 

300 

56.75 

24.25 

21 

300 

58.30 

25.00 

2 

373 

59 .  .35 

20.38 

21 

523 

65.17 

15.96 

21 

609 

67.73 

14.49 

21 

II.   THEORETICAL  CONSIDERATIONS 

A.   CRYSTALLINE  EFFECTS  ON  NEUTRON  SCATTERING 

Since  the  wavelengths  of  thermal  neutrons  are  of  the  same  order  of 
magnitude  (about  1.8  A  at  room  temperature)  as  the  interatomic  plane 
spacing  for  crystalline  materials  such  as  graphite,  certain  optical 
effects  are  observed.   That  is,  all  the  interference  effects  associated 
with  light  and  x-rays  can  be  duplicated  by  the  neutrons.   It  will  be 
shown  that  there  are  energy  regions  where  these  optical  effects  may  or 
may  not  be  present.   When  they  are  present  we  have  coherent  scattering; 
that  is,  the  kind  of  scattering  that  produces  a  wave  capable  of  inter- 
fering with  the  incident  neutron  wave.   In  the  absence  of  optical 
effects,  the  scattered  wave  does  not  interfere  with  the  incident  wave, 
resulting  in  "diffuse"  scattering,  called  incoherent  scattering. 

The  coherent  scattering  may  be  composed  of  both  elastic  and .. 
inelastic  scattering.   Incoherent  scattering  may  also  be  either  elastic 
or  inelastic,  broken  down  generally  as  follows:   elastic  incoherent  may 
be  due  to  the  presence  of  different  isotopes  in  the  atomic  lattice 
(which  gives  diffuse  scattering) ,  or  due  to  the  neutron  spin  relative 
to  the  nuclear  spin;  inelastic  incoherent  which  is  due  to  the  energy 
exchange  between  a  neutron  and  the  medium's  lattice  vibrations,  with 
the  neutron  absorbing  one  or  more  phonons  (the  lattice  vibration 
quantum)  of  energy.   This  phonon  energy  may  be  from  vibrations  in  i.    y 
of  the  degrees  of  freedom. 

The  total  cross  section  may  then  be  written  for  the  thermal  region 

as: 

0   =  o  +  a.   +  o     +  o , . 
t      a    1    o    dis 


where  o   is  the  capture  cross  section,   a   is  the  coherent  elastic 
a  r  o 

scattering  cross  section,  a,.        is  the  incoherent  elastic  scattering 
cross  section  due  to  isotope  presence  or  neutron/nucleii  spin,  and 
o1   is  the  inelastic  incoherent  and  coherent  scattering  cross  section 
due  to  phonon  interchanges.   In  the  thermal  region,  only  the  one  phonon 
effect  (neglecting  the  two  or  more  phonon  effects)  makes  any  contribu- 
tion to  inelastic  scattering  113]. 

In  the  case  of  graphite  some  of  these  terms  may  be  eliminated. 
Since  coherent  inelastic  scattering  is  low  compared  to  the  coherent 
elastic  scattering  it  is  generally  neglected  at  energies  below  1.0  ev. 
The  capture  cross  section  of  graphite  is  small  but  proportional  to  1/v. 
Since  graphite  is  generally  monoisotopic  and  has  zero  spin,  the 
incoherent  elastic  scattering  term  drops  out.   Therefore  there  exists 
only 

cr   =  a  +  a  (above  Bragg  cut-off) 

t       o     a  &e>         / 


or 


a       =  a.  +  a  (below  Bragg  cut-off) 

£  X        3. 


In  coherent  scattering  the  mechanism  by  which  waves  capable  of 
interfering  with  the  incident  neutron  wave  is  governed  by  Bragg' s  law, 
written  as 

nX  =  2d  sin  6 

where  n  is  the  order  of  reflection,  X  is  the  wavelength  corresponding 
to  the  neutron  energy,  d  is  the  plane  spacing  and  0  is  the  angle 
between  the  direction  of  neutron  propogation  and  the  crystal  plane. 
From  this  basic  equation  it  is  seen  that  there  are  many  wavelengths 
(or  energies)  for  which  the  coherent  scattering  mechanism  is  possible, 


based  on  many  different  plane  spacings.   There  is  also  a  minimum  energy 
corresponding  to  a  maximum  wavelength,   X  A  =  2d  ,  below  which 
coherent  scattering  is  not  possible  for  any  given  plane  spacing.   For 
the  largest  plane  spacing,  the  002  plane  in  graphite,  this  minimum 
energy  is  called  the  Bragg  cut-off,  and  no  further  coherent  scattering 
can  take  place  at  lower  energies.   For  graphite,  the  average  room 
temperature  crystalline  002  plane  spacing  d  is  3.3538  A  (Angstroms). 
Using  the  deBoglie  relation 

h       h 


mv 


h 

v  "  Am" 


where  h  is  Planck's  constant,  p  is  momentum,  m  is  the  neutron  mass  and 
v  is  the  neutron  velocity,  the  energy  of  this  Bragg  cut-off  can  be 
determined. 


KE  =   y  mv 


ra  ■  -b<!>2 


KEMIN      2m  (  2d  } 


Substitution  of  the  values  for  neutron  mass,  Planck's  constant  and 
plane  spacing  yields  a  Bragg  cut-off  energy  of  0.001823  ev.   Thus, 
for  all  energies  above  this  value  (a  region  occupied  by  neutrons  of 
the  thermal  group)  coherent  elastic  scattering  is  the  dominant  factor, 
while  at  lower  energies,  for  neutrons  of  the  cold  group,  incoherent 
inelastic  scattering  alone  dominates. 

Although  the  frequency  of  the  lattice  vibrations  is  largely 
independent  of  temperature,  the  amplitude  increases  with  increasing 


temperature  [8],   This  increased  amplitude  results  in  increased  lattice 
energy  and  therefore  larger  one  phonon  energies,  giving  greater 
incoherent  inelastic  scattering.   Additionally  it  can  be  shown  [10] 
that  the  incoherent  inelastic  scattering  is  proportional  to  the  wave- 
length, or  follows  the  "1/v  law",  and  that  it  increases  rapidly  with 
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temperature  (about  linear  by  observation,  as  opposed  to  the  T  behavior 

predicted  by  the  usual  Debye  frequency  spectrum) .   Capture  cross 
section,  however,  is  temperature  independent  for  all  thermal  energies. 
Figures  1  and  2  show  the  coherent  and  incoherent  scattering  cross 
section  regions.   The  combined  effects  of  these  graphite  scattering 
cross  sections,  along  with  its  capture  cross  section  is  seen  in 
Figure  3. 

B.   DIFFUSION  OF  NEUTRONS  IN  A  CRYSTALLINE  MODERATOR 
The  diffusion  length  can  be  defined  by 


T2     1   ,  .     1 

L    =  3  •  \x  •  r 

a 


L2     * 


3   ncr     I 

tr    a 


where  A    is  the  transport  mean  free  path,   Z   is  the  macroscopic 

capture  cross  section,  n  is  the  number  of  atoms  per  unit  volume  and 

o   is  the  transport  cross  section.   In  the  thermal  neutron  group 

region,  above  the  Bragg  cut-off,  the  value  of  a         will  be  largely 

that  of  the  elastic  coherent  scattering  cross  section,  o   .   It  is 

o 


generally  assumed  to  be  given,  away  from  the  immediate  vicinity  of  the 

2 
rs  ^"  "  3A*  *1-',  where  °s 


2 
Bragg  cut-off,  by  its  asymptotic  value  a   (1  -  -r-r-)  11],  where  a 
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Elastic  (Bragg)  coherent  scattering  cross  section  of  graphite,  [6] 
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Inelastic  incoherent  scattering  cross  section 
of  graphite  (k  is  Bolzmann's  constant),  [1] 
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is  the  microscopic  (coherent  elastic)  scattering  cross  section  and  A 
is  the  atomic  mass  number.  The  temperature  dependence  of  this  cross 
section  is  very  small.   In  the  cold  neutron  group  region  a  can  be 

taken  to  be  equal  to  o.  ,  the  incoherent  inelastic  cross  section  due 
to  the  one  phonon  effect.   As  pointed  out  above,  and  shown  in  Figure  2, 
this  cross  section  increases  with  temperature.   Therefore,  the  diffusion 
length  should  decrease  with  increase  of  temperature  for  energies  in  the 
cold  neutron  group,  while  remaining  essentially  constant  with  tempera- 
ture for  energies  above  the  Bragg  cut-off. 

The  time  independent,  isotropic  diffusion  equation  for  thermal  flux 
can  be  written  as 


V2Hr)  -  -^  4>(r)  =  0  (1) 

V 


where 


D    1/2 
K     =   (  -=  ) 

a 


and  where  D  and  E  are  the  average  values  in  the  thermal  region  for 
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diffusion  coefficient  and  macroscopic  capture  cross  section,  respectively. 
Equation  (1)  can  be  written  for  two  regions  of  energy,  the  cold  and 
thermal  groups,  with  the  values  of  D  and  E   being  determined  over 

3. 


the  appropriate  energy  range. 


,ET~lev 


Dl  = 


f 

J      D(E)  <},(E)  dE 
0.0018 


ET~lev 


f 

J  HE)  dE 

0.0018 
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Thus, 


and 


/  0.0018 
D(E)  ME)  dE 


D2       0.0018 


J         ME)  dE 

•'o 


I 


ET~lev 


I    (E)  ME)  dE 
0.0018   a 

al  =    v     i 

i«ET~lev 

I         ME)  dE 

•^0.0018 


I 


0.0018 


I    (E)  ME)  dE 

Ea.  -  -5 : 

a2 

J  .0.0018 
♦(E)  dE 
0 


vY(r)  -  -^-2  Mr)   -  0  Cl-a) 

LT1 


V2^>2Cr)  -  -^-2  $2(r)   =  0  (1-b) 

LT2 


where  the  subscripts  1  and  2  refer  to  the  thermal  neutron  group  and 
cold  neutron  group,  respectively.  With  the  source  placed  at  the  center 
of  the  base  of  a  rectangular  parallelopiped,   <J)(x,y,0)   will  be 
symmetric  about  the  origin  and  the  flux  will  be  an  even  function  of 


14 


both  x  and  y  throughout  the  pile.   The  solution  to  equations  (1-a)  and 
(1-b)  has  the  same  form.   It  can  be  shown  that  (dropping  subscripts  for 
simplicity)  the  solution  is  I15J 


.  /     n      t»  a     •  i-     /   \       mirx       nuy    /ON 

(bCx.y.z)   =   E    A    smh  y   Cc-z)  •  cos  •  cos  — *-        (2) 

r  mn       mn            a        a 
m,n 

odd 


where 

a/2  ^3/2 


mn 


n 

J-a/2j-a 


miTx    nu 


1 


A    =  — x g     I      4>_(x,y,0)cos cos—- x  dx  dy 

mn      2   .  ,.  .f     |      X   "       a      a      y 

a  sxnh(Y   c) 


-a/2*J-a/2 


Y   is  the  relaxation  length  and  is  expressed  as 
mn  ° 


Y   2  o  _J_   (S]L)2    .nrr.2  (3) 

Ymn      .  2  +  ^  a;    V  a; 
LT 


c  is  the  vertical  extrapolated  height  of  the  pile. 
a  is  the  extrapolated  horizontal  dimension  of  the  pile. 
x,y,z  are  the  co-ordinates  from  an  origin  at  the  bottom  center  of 
the  pile. 
Since  measurements  are  made  along  the  centerline  of  the  pile 
where  x  =  y  =  0  ,  equation  (2)  simplifies  further  to 


<b(0,0.z)      I       A   sxnh(c-z) 
T  mn 

m,n 

odd 

or  in  exponential  terms 

,  Y   Cc-z)    ~y   (c-z) 

♦  CO.O.z)   =  \       I         Amn  I  e  ffln      -  e  mn      ]        CA) 
m,n 
odd 
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In  a  pile  measurement  the  height  c  is  much  larger  than  the  diffusion 

length  so  that  in  areas  not  too  near  the  top,  the  term  Y   (c-z) 

inn 

becomes  very  large  and  the  second  term  of  equation  (4)  can  be  neglected, 

Then 

n  y   c    -y   z 

$(0,0, z)   =  tt       I       A   e     *  e 
2       mn 
m,n 

odd 


and  by  defining 


C    =  r?  A   e 

mn     2  m,n   mn 

odd 


-Y   z  (5) 


d>(0,0,z)  =    I       C   e  mn 

mn 
m,n 

odd 


~Y   z        -y   z        -y   z 
<})(0,0,z)   =  C  -  e      +  C^  e      +  C   e      +  .  .  . 

Hughes  [10]  points  out  that  as  the  distance  from  the  source  increases, 
the  effect  of  the  higher  harmonics  (Y-io>  Yot  >  et  cetera)  becomes 
negligible  and  that  the  neutron  distribution  is  practically  exponential 
beyond  that  distance.   Also  Kaiser  and  Kimel  [12]  have  computed  the 
relative  contributions  of  the  y-io  harmonic  for  a  graphite  pile  as  a 
function  of  distance  from  the  source.   They  show  Y-io   to  ^e  approxi- 
mately one-tenth  the  contribution  of  Y-i-i   f°r  a  constant  thermal 
source  at  a  distance  of  20  inches  from  the  source. 

Therefore  by  concentrating  on  the  mid-axis  range  —  from  35  cm  to 
80  cm  —  it  is  possible  to  neglect  the  higher  harmonics  and  end  effects 
so  that 
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or 


«"Y   Z 

*(0,0,z)  ;=£  Cu  e  U 


In  $  ■  In  Cn  -  y-^z  C6) 


which  is  a  simple,  one  degree  polynomial  equation  of  intercept   In  C. 
and  slope  -y,.   on  a   In  <j>  versus   z  plot.   The  form  of  this 
equation  is  valid  for  both  thermal  and  cold  neutron  groups. 

Least  squares  polynomial  approximations  over  the  region  of  validity 
for  each  group  results  in  the  determination  of  their  respective  Y-ii* 
Since  the  relaxation  length  is  a  function  of  L  and  it  has  been  shown 
that  diffusion  length  for  an  infinite  medium  changes  across  the  Bragg 
cut-off,  it  follows  that  the  values  of  Y-i-i*  or  the  slopes  of  the  plot, 
should  change  also.   This  will  occur  at  some  point  in  the  pile  at  which 
the  cold  neutrons  have  been  filtered  through  in  appreciable  number. 

Similarly  the  degree  of  slope  in  each  instance  will  change  with  a 
temperature  change.   The  greater  degree  of  change  will  be  present  in 
the  cold  group  region  neutrons,  that  is,  the  region  where  the  number 
of  cold  neutrons  is  appreciable.   Since,  as  was  pointed  out,  the  inco- 
herent inelastic  scattering  in  this  region  is  highly  temperature 
dependent,  an  increase  in  temperature  will  result  in  increased 
scattering  cross  section  and  reduced  diffusion  length.  A  lower 
diffusion  length  is  associated  with  an  increased  slope  Yii   (refer  to 
equation  (3)).   In  the  thermal  neutron  group,  the  transport  mean  free 
path  has  been  found  by  Lloyd,  et  al.  [16]  to  decrease  slightly  with 
temperature.   Therefore  with  higher  temperatures  the  diffusion  length 
would  increase  and  the  slope  Y-,-,  would  decrease.   The  overall  effect 
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of  high  temperature  will  then  be  to  increase  the  slope  of  the  first 
part  of  the  neutron  distribution  curve  (before  the  "knuckle")  and 
decrease  the  slope  of  the  second  (or  cold  neutron)  part  of  the  curve 
(beyond  the  "knuckle") .   At  a  sufficiently  high  temperature  the 
"knuckle"  will  vanish  and  a  single  slope  will  exist  for  the  entire 
thermal  range.   According  to  cross  section  data  from  Ref.  {11],  this 
temperature  should  be  higher  than  1020°K. 
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III.   PROJECT  DESCRIPTION 

It  was  proposed  to  verify  experimentally:   first,  that  the  thermal 
neutron  diffusion  lengths  above  and  below  the  Bragg  cut-off  energy  did, 
in  fact,  differ;  and  second,  that  this  difference  decreased  as  tempera- 
ture increased.   The  neutron  flux  plot  should  change  from  a  distinct 
"knuckle"  at  the  location  where  the  number  of  cold  neutrons  becomes 
appreciable  to  a  single,  continuous  line  at  very  high  temperatures. 
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IV.   TEST  APPARATUS 

A.   DESCRIPTION 

1.  General 

It  was  originally  planned  to  utilize  the  AGN  201  reactor  at  the 
Naval  Postgraduate  School  reactor  facility  as  the  neutron  source  for 
this  experiment.   However,  the  physical  construction  of  the  reactor  did 
not  provide  for  removal  of  any  segments  of  its  shielding.   It  would 
have  been  necessary  to  construct  the  graphite  pile  outside  the  reactor  — 
a  distance  of  some  7.5  feet  from  the  reactor  core  center  —  over  one  of 
the  access  holes  or  the  "glory  hole".   It  would  be  necessary  for  the 
reactor  to  be  operated  at  its  maximum  power  of  1000  watts  (which  could 
only  be  sustained  for  about  ten  minutes)  in  order  for  sufficient  neutron 
flux  to  be  obtained  in  this  manner.   This  short  time  period  was  quite 
insufficient  for  any  flux  measurement  in  the  proposed  sigma  pile. 

Therefore  it  was  decided  to  utilize  three  Plutonium-Beryllium 
neutron  sources:   two  sources  of  1.8  x  10  n/sec  each,  and  one  source 
of  5.4  x  10  n/sec.   These  were  the  only  other  continuous  neutron 
sources  available.   The  most  probable  energy  of  these  neutron  sources 
is  about  4.0  Mev  [3],  so  it  would  be  necessary  to  provide  for  additional 
slowing  down  time. 

2.  Graphite  Pi] a 

The  graphite  pile  consisted  of  a  rectangular  parallelopiped 
constructed  of  one  hundred  twenty  eight  (128)  graphite  bricks,  each  one 
3  x  6  x  12  inches.   The  base  was  two  feet  square  and  the  height  was 
four  feet.   Figure  4  gives  a  cross  sectional  view  of  the  pile's  general 

arrangement.   Great  Lakes  Carbon  Corporation,  Grade  HC,  graphite  of 

3 

density  1.60  gm/cm  was  utilized.   The  use  of  this  relatively  impure 
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graphite  as  opposed  to  AGOT  nuclear  graphite  was  considered  acceptable 
in  view  of  the  experiment's  primary  aim  in  investigating  temperature 
dependence  rather  than  absolute  values.   A  three  sixteenth  inch  vertical 
hole  was  drilled  down  the  entire  four  foot  centerline  of  the  pile  to 
provide  for  the  detecting  wire  used. 
3.   Heating  and  Insulation 

Heating  of  this  pile  was  achieved  by  twenty  calrod  strip  heaters 
of  21.5  inches  heating  length,  1.5  inches  width  and  3/8  inch  thickness 
arranged  about  all  four  sides,  the  top  and  bottom.   Close  fit  grooves, 
1.5  by  3/8  inch,  were  cut  in  the  graphite  bricks,  three  inches  in  from 
the  side  faces  and  six  inches  in  from  the  top  and  bottom,  prior  to 
assembly  to  accomodate  these  heaters.   The  grooves  were  a  full  24 
inches  long,  the  same  size  as  the  pile,  in  order  to  provide  for  thermal 
expansion  of  the  heater  elements.   Figure  5  shows  the  external  arrange- 
ment of  graphite  bricks  and  groove  openings.   Each  heater  was  rated  for 
750  watts  at  220  volts,  providing  for  15  kVA  of  heating  capacity. 
Heaters  were  regulated  in  pairs  by  ten  variacs  which  provided  excellent 
control  of  heat  throughout  the  pile.   Figure  6  shows  control  table 
arrangement  with  the  variacs  and  thermocouple  readout  instruments. 

The  design  of  insulation  was  entirely  dictated  by  the  availa- 
bility of  insulating  materials  within  the  campus.  Fire  brick  was  used 
as  the  base  insulation  since  it  was  strong  enough  in  compression  to 
support  the  1590  pound  weight  of  the  graphite  and  also  had  excellent 
insulating  properties.   A  sheet  of  three  sixteenth  inch  asbestos  cement 
board  separated  the  fire  brick  from  the  half  inch  steel  plate  which 
acted  as  a  table  to  support  the  pile  over  the  neutron  source.   Both  the 
steel  plate  and  fire  brick  had  a  three  inch  hole  at  the  pile  centerline 
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Fig.  6 
Control  Table  Arrangement 

to  reduce  neutron  scattering  in  that  immediate  vicinity.   Two  types  of 
insulating  block,  using  readily  available  material,  were  installed 
about  the  four  sides  at  a  distance  of  about  three  inches  from  the  pile. 
A  two  inch  layer  of  glass  wool  insulation  attached  to  the  underside  of 
a  plywood  board  encased  the  top  of  the  insulating  walls,  also  with  an 
air  gap,  thus  completely  insulating  the  graphite  pile.   The  air  gaps 
were  provided  in  an  attempt  to  reduce  to  a  minimum  the  backscatter  that 
may  have  occurred  from  the  insulating  material,  thus  making  possible 
the  use  of  the  simplified  diffusion  equation  given  as  equation  (1) . 
Figure  7  shows  the  lower  corner  of  the  graphite  pile  to  illustrate  the 
base  and  side  insulation,  as  well  as  the  heating  element  connections. 
A.   Temperature  Determination 

Ten  copper-cons tantan  miniature  thermocouples  were  arranged  in 
groups  of  five  each  in  planes  that  were  twelve  and  thirty  six  inches 
respectively  from  the  base.   In  each  plane  one  thermocouple  was  at  the 
centerline  of  the  pile  and  the  other  four  were  three  inches  from  the 
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centerline.   The  closest  any  thermocouple  was  to  any  heating  element 
was  6.72  inches  and  the  most  distant  was  7.51  inches.  Due  to  symmetry 
of  all  heating  element  locations,  the  thermocouples  were  equally 
affected  by  the  aggregate  heaters.  This  arrangement  of  thermocouples 
allowed  determination  of  temperature  uniformity  throughout  the  volume 
of  graphite  for  which  the  diffusion  length  was  measured.   A  multiple 
switch  was  utilized  to  give  a  common  reference  junction  of  0°C  to  all 
thermocouple  readings.   A  plain  water  ice  bath  in  a  thermos  jug  served 
as  the  reference  junction.  A  digital  potentiometer  provided  the  voltage 
readout  continuously  to  within  a  tenth  of  a  millivolt.   A  galvanometer- 
type  millivolt  potentiometer  with  accuracy  to  31.0  microvolt  was 
utilized  for  final  voltage  determination  after  the  system  had  settled 
out  and  for  accurate  temperature  determination  during  irradiation  runs. 
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It  was  recognized  that  copper  was  probably  the  worst  possible 
constituent  to  use  in  a  nuclear  radiation  experiment  for  temperature 
determination.   As  pointed  out  by  Ross  118] ,  the  percentage  change  of 
resistivity  with  temperature  in  a  neutron  flux  can  be  expressed  by: 


%  &&■     =   220  x  10"22  At 
po 


compared  to  that  for  constantan: 

%  &■     =  4.7  x  10"22  4>t 
Po 

where  Ap   is  the  change  of  resistivity,  p   is  the  original  resistiv- 
ity,  <f>  is  the  neutron  flux  and  t  is  the  time  in  the  neutron  flux  in 
seconds.  However,  the  low  neutron  flux  available  for  this  experiment 
and  the  relatively  short  exposure  times  permitted  the  neglect  of  this 
shortcoming.   As  with  the  decision  to  use  plain  water  versus  pure  water 
in  the  ice  bath,  the  accuracy  of  termperature  determination  was  deemed 
sufficient  to  within  five  degrees  centigrade. 
5.   Neutron  Source 

The  three  Plutonium-Beryllium  (Pu-Be)  sources  were  grouped 
together  in  a  horizontal  plane  beneath  the  base  of  the  pile  directly 
under  the  centerline  axis.   A  single  one  and  three  quarter  inch  layer 
of  paraffin  was  placed  above  the  sources  to  aid  in  slowing  down  the 
high  energy  neutrons  before  they  entered  the  graphite.  Numerous 
arrangements  of  source  height  position  and  thickness  of  paraffin  were 
tried  before  the  above  set  up  was  decided  upon.   It  was  found  that  the 
paraffin  not  only  slowed  down  but  absorbed  the  neutrons  so  much  that 
the  flux  was  reduced  to  an  unacceptable  level.   Therefore  only  one 
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layer  was  utilized.   This  layer  and  the  source  under  it  were  placed  as 
close  to  the  base  as  possible  since  it  was  also  postulated  that  the  air 
gaps  present,  and  the  different  medium  through  which  the  neutrons 
travelled  before  reaching  the  graphite  all  contributed  to  the  low  flux 
within  the  graphite.   Since  these  neutrons  were  of  a  higher  energy  than 
normally  encountered  in  experiments  of  this  nature,  it  was  thought  that 
they  would  require  a  much  larger  amount  of  moderating  material  to  reach 
thermal  energies.   However  a  quick  calculation  of  the  number  of  collisions 
required  to  achieve  this  energy  in  hydrogenous  moderators  (assuming  the 
same  logarithmic  decrement  as  water)  showed  as  follows : 


/  n  m      .  1  ,   Eo  .    1  ,   4  x  106 

4.0  Mev  neutrons:    n  =  —  In  —  -  ~  no  In 


E     0.92       1 
n  =  16.5  collisions  to  reach  1.0  ev 

oa«                   l-.o_l-.2x  10 
2.0  Mev  neutrons:    n  =  —  In  —  -  A  no  In  z 

n  =  15.78  collisions  to  reach  1.0  ev 

With  a  transport  mean  free  path  of  0.395  cm  for  paraffin  [20]  (the 
average  distance  between  collisions) ,  this  indicated  rhnt   only  an 
additional  0.28A  cm  would  be  required  to  achieve  thermal  energies. 
Blocks  of  paraffin  surrounded  the  sides  of  the  source  to  a 
thickness  of  two  and  one  half  feet  and  were  arranged  outside  the 
graphite  insulation  to  protect  personnel  in  the  area.   The  neutron 
leakage  through  this  arrangement  was  surveyed  and  found  negligible. 
6.   Neutron  Detection 

It  was  originally  planned  that  a  BF_  neutron  detector,  in 
conjunction  with  a  multi-channel  analyzer,  would  be  used  to  take  an 
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integrated  count  at  each  axial  position.   Investigation  showed  that  this 
detector  would  withstand  temperatures  only  as  high  as  150°C.   Similar 
restrictions  were  found  to  exist  for  solid  state  detectors  as  well  as 
indium  and  cadmium,  which  melt  at  156°C  and  320. 9 °C,  respectively.   A 
few  solid  state  detectors  were  available  commercially  but  to  achieve  the 
sensitivity  required  for  the  low  flux  measurements  in  this  experiment, 
their  size  was  much  too  large.   Gold  foils,  though  commonly  used  for 
high  temperature  measurements  offered  such  a  long  half-life  as  to  be  not 
feasible  from  irradiating  time  considerations.   A  detecting  wire  of  11.4 
per  cent  dysprosium-aluminum  alloy  was  finally  chosen,  based  on  its  high 
melting  temperature  and  convenient  (2.35  hour)  half-life.  A  purity  of 
99.94  per  cent  was  certified  by  the  manufacturer  and  impurities  listed 
by  name  indicated  no  need  for  their  consideration  or  possible  effect  on 
the  dysprosium  capture  cross  section;  although  this  factor  was  taken 
into  account  in  computing  the  total  number  of  dysprosium  atoms,  described 
in  Appendix  A.   An  average  atomic  weight  of  162.5  was  determined  based 
on  the  percentage  of  the  various  isotopes  to  be  found  in  natural 
dysprosium.   A  value  of  930.0  barns  was  used  as  the  capture  cross 
section,  as  recommended  [7].   This  coincides  with  the  value  obtained  by 
weighting  the  cross  sections  of  each  isotope  in  accordance  with  their 
abundance  in  natural  Dy,   Dysprosium  has  a  small  response  to  epithermal 
neutrons  [14,19],   The  maximum  epithermal  contribution  to  the  total 
dysprosium  activity  has  been  found  to  be  less  than  two  per  cent  [4J. 
Therefore  the  activity  determined  with  this  dysprosium  wire  was  consid- 
ered to  be  completely  due  to  thermal  flux.   The  detecting  wire  was 
positioned  in  the  three  sixteenth  inch  axial  hole  along  the  entire  pile 
height  during  irradiation.  A  thin  glass  tube  held  the  wire  segments  in 
position. 
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On  initial  irradiation  tests  the  reactor  facility  counting 
equipment  was  utilized.   This  consisted  of  a  thallium-activated  sodium 
iodide  crystal  and  photo-multiplier,  a  pulse  height  analyzer  for  total 
decay  count  determination,  and  a  multi-channel  analyzer  for  energy 
spectrum  readout  and  analysis  (as  a  side  benefit  and  not  directly  con- 
nected with  the  object  of  the  experiment).   It  was  later  found  to  be 
more  practical  to  utilize  the  counting  equipment  of  the  laboratory  adja- 
cent to  the  experimental  room.   This  involved  a  Geiger-Mueller  tube 
operating  with  1400  volt  D.C.  bias,  and  pulse  height  analyzer.   The 
efficiency  of  the  G.M.  counting  system  for  beta  counting  was  determined 
to  be  approximately  twenty-seven  per  cent,  while  that  of  the  photo- 
crystal  for  gamma  counting  in  the  reactor  facility  was  about  twenty-nine 
per  cent,  so  that  no  appreciable  advantage  could  be  gained  by  the 
latter' s  utilization  other  than  familiarization  with  another  piece  of 
equipment. 

B .   OPERATION 

Test  runs  were  made  at  room  temperature  and  at  an  elevated  tempera- 
ture of  481°K.   The  following  operational  procedures  were  followed  in 
both  cases  except  for  temperature  control. 

Two  main  fifty  ampere  circuit  breakers  were  turned  on  to  provide  a 
maximum  of  6.8  amperes  to  each  of  ten  variacs,  which  in  turn  provided 
3. A  amperes  to  each  of  two  heating  elements.   The  variacs  were  connected 
such  that  five  of  them  controlled  the  bottom  ten  heaters  and  five  of 
them  controlled  the  top  ten  heaters.   Similarly  each  variac  controlled 
the  same  two  elements  at  a  particular  height  so  that  even  control  of 
temperature  could  be  achieved.   All  variacs  were  gradually  increased  to 
slowly  raise  the  pile  temperature  to  the  desired  level  with  the  minimum 
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amount  of  "hunting".   It  was  found  that  the  temperatures  at  each  thermo- 
couple grouping  were  consistently  even,  within  2°C  of  each  other.   It 
was  thus  only  necessary  to  adjust  the  top  with  the  bottom  groupings. 
This  was  a  fairly  easy  task  and  temperature  variations  between  any  of 
the  ten  thermocouples  amounted  to  a  maximum  of  plus  or  minus  4.56°C. 
The  final  test  temperature  value  is  the  mean  of  the  ten  readings,  which 
are  about  constant,  over  the  entire  irradiation  time  span. 

Once  the  specified  temperature  equilibrium  was  achieved,  the  three 
sources  were  placed  in  position  under  the  graphite  pile.   In  order  to 
utilize  the  counting  equipment  and  obtain  the  activity  at  a  particular 
location  it  was  necessary  to  cut  the  detecting  wire  into  one  inch 
segments.   Due  to  the  cost  of  the  wire,  segments  had  to  be  used  over 
again,  necessitating  a  requirement  for  being  able  to  insert  and  remove 
segments  to  a  known  axial  position.   This  was  achieved  through  the  use 
of  a  thin  glass  tubing  the  length  of  the  graphite  pile,  into  which  the 
wire  segments  could  be  inserted.   The  bottom  of  the  tube  was  sealed  with 
a  glue  to  prevent  segments  from  dropping  out.   A  maximum  of  eight  seg- 
ments could  safely  be  inserted  in  the  glass  tubing  without  the  glass 
fracturing.   Since  48  data  points  were  being  taken  this  meant  at  least 
six  irradiation  periods  for  each  test  run.  An  average  irradiation 
period  of  eight  hours  provided  ninety  per  cent  of  full  saturation. 
After  irradiation  the  tube  was  withdrawn,  the  bottom  two  elements 
removed  by  snapping  off  the  glass  tubing  at  that  point,  adding  a  drop 
more  glue  to  the  new  bottom  of  the  tube,  and  reinserting  the  tube, 
maintaining  segment  axial  position,  to  continue  the  irradiation  proce- 
dure while  the  two  elements  removed  were  being  counted  for  decay 
activity,   Due  to  the  low  flux,  long  (thirty  minute)  counting  times  per 
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segment  were  required.   Since  even  during  the  short  time  of  absence  from 
the  flux  for  segment  removal  a  significant  decay  occurred  in  the  re- 
maining segments,  it  was  necessary  to  derive  the  iterative  formula  for 
the  activity  at  each  subsequent  removal  time  as  described  in  Appendix  A. 

Upon  removal  of  the  two  segments  a  period  of  seven  minutes  was 
allowed  to  elapse  in  order  to  reduce  the  75  second  half-life  of  Dy 
to  a  negligible  effect.   During  this  time,  ten  minute  background  counts 
were  obtained,  to  be  averaged  together  and  applied  to  the  segment  total 
activity  count.   During  the  initial  stages  of  the  experiment,  using  the 
reactor  facility  equipment,  the  energy  spectrum  of  the  decaying  wire 
segment  was  determined  from  the  multi-channel  analyzer.   Although  the 
decay  rates  were  very  small,  a  typical  dysprosium  gamma  ray  spectrum 
was  obtained  with  no  indication  of  any  other  material  being  present, 
validating  the  stated  composition  of  dysprosium  and  aluminum.   Through- 
out the  entire  experimental  testing,  accurate  determination  and  recording 
of  time  was  maintained. 

The  individual  wire  segments  were  used  in  the  same  axial  position 
for  both  runs  to  insure  consistency  in  obtaining  data. 

C.   FAILURES  AND  MALFUNCTIONS 

The  scope  of  the  original  project  definition  had  to  be  curtailed 
drastically  for  a  number  of  reasons,  some  already  stated: 

1.   Detector  Restrictions 

The  high  temperatures  sought  to  be  achieved  drastically 
restricted  the  choice  of  detectors  available.   The  state  of  the  art  in 
neutron  detectors,  especially  in  high  temperature  ranges,  is  geared 
toward  very  high  flux  ranges  as  found  in  reactors.   This  required  the 
use  of  a  dysprosium  wire  which,  while  offering  negligible  self-shadowing 
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and  being  sensitive  mostly  to  thermal  neutrons,  also  was  too  small  to 
make  up  the  adverse  effects  of  low  flux.   Since  the  wire  was  very 
expensive  it  was  necessary  to  find  some  method  of  re-using  the  segments 
after  they  had  been  previously  used. 

A  tubing  arrangement  was  decided  upon  as.  offering  the  simplest 
method  of  positioning  and  handling  the  segments.   Glass  tubing  was 
readily  available  but,  due  to  its  brittleness,  there  was  a  high  proba- 
bility of  tube  fracture  and  subsequent  loss  inside  the  pile.   This 
method  was  to  be  used  only  as  a  last  resort.   Three  different  metal 
tubes  of  various  sizes  were  tried.   However,  in  each  case  the  segments 
bound  up  inside  the  tubing,  necessitating  the  cutting  of  the  tubes  to 
remove  the  segments.   It  was  finally  resolved  that  the  wire  was  too 
soft  to  be  stacked  inside  a  tube  or  pushed  through  it  by  a  rod.   The 
slightest  buckle  produced  a  jamming  against  the  tube  walls,  inducing 
further  binding.   Therefore  the  glass  tubing  was  resorted  to,  as 
described  in  the  operating  procedures,  with  utmost  caution  being  taken 
to  avoid  fracture  inside  the  pile. 
2.   Heat  Considerations 

The  Pu-Be  source  provides  neutrons  from  the  alpha  particles 
(given  off  by  plutonium  decay)  interacting  with  the  beryllium  to  give 
off  neutrons  (Be(a,n)  reaction).   There  is  the  possibility  that  some 
alpha  particles  would  pick  up  electrons  to  become  helium  gas  before  an 
interaction  with  beryllium  could  occur.   Thus  it  was  decided  to  keep 
the  sources  at  room  temperature  to  prevent  a  possible  rupture  of  the 
source  cladding  due  to  gaseous  expansion,  despite  assurances  from  the 
manufacturer  that  no  hazard  existed.   Coincidentally,  by  maintaining  the 
sources  at  one  temperature  for  all  tests  a  more  realistic  comparison 
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of  data  would  be  possible  since  the  source  spectrum  could  be  a  function 
of  temperature. 

This  tremendously  complicated  the  construction  of  the  experi- 
mental set-up.   With  the  existing  insulation  installed  below  the  base, 
a  design  was  provided  to  have  an  air  channel  over  the  source  and  its 
moderating  paraffin  layer.   A  small  blower  motor  was  to  have  been 
installed  to  force  cooling  air  through  the  channel  and  out  an  exhaust 
ducting.   However,  when  it  became  apparent  that  the  flux  level  was 
going  to  be  too  low,  it  was  necessary  to  do  away  with  the  air  channel 
and  move  the  paraffin  and  source  as  close  to  the  base  as  possible.   A 
remote  reading  bourdon  temperature  gage  was  installed  between  paraffin 
layers  in  order  to  keep  the  source  temperature  under  constant  surveil- 
lance. 

During  heating  operations  the  temperature  variation  between  all 
measuring  thermocouples  was  easily  kept  within  the  limits  desired.   It 
was  found  however  that  over  a  period  of  time  the  steel  plate  began  to 
heat  up  to  such  a  degree  as  to  affect  the  paraffin.   This  was  traced 
to  the  three  inch  square  hole  in  the  firebrick  below  the  graphite. 
Although  there  was  also  a  hole  in  the  steel  plate  at  that  location, 
the  single  sheet  of  asbestos  sheet  cement  did  not  offer  long  term  heat 
insulation  properties  and  conducted  the  heat  throughout  the  steel  plate 
underneath.   By  filling  this  hole  in  the  firebrick  arrangement  with 
plastic  refractory  and  a  small  sheet  of  insulating  block,  it  was 
possible  to  eliminate  the  problem  of  the  hot  steel  plate.  As  an 
additional  measure,  some  of  the  surrounding  paraffin  in  close  proximity 
to  the  steel  plate  was  removed,  while  leaving  the  paraffin  around  the 
source.   This  provided  a  cavern  effect  around  the  center  column  of 
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paraffin  and  source  (which  remained  in  place,  close  to  the  graphite 
base) .   The  air  blower  originally  planned  for  was  then  installed  to 
circulate  air  through  this  cavern  and  around  the  source,  rather  than 
over  it  as  was  originally  planned.   This  proved  quite  satisfactory  and 
no  further  heat  problems  were  encountered. 
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V.   DATA  ANALYSIS  AND  RESULTS 

A.   DATA  REDUCTION 

A  total  of  forty  eight  data  points  were  collected  for  the  room 
temperature  test.   After  the  break  in  the  slope  was  located  from  this 
run,  only  the  center  twenty  two  positions  were  investigated  for  the 
high  temperature  test.   The  saturation  activity  and  corresponding 
thermal  flux  was  computed  as  shown  in  Appendix  A.   The  uncertainties 
are  described  in  Appendix  B. 

The  natural  log  of  the  flux  was  then  plotted  versus  the  height  in 
order  to  determine  where  the  break  in  the  slope  might  occur.  As  noted 
in  the  theory,  only  the  center  portion  of  the  pile  was  to  be  investi- 
gated.  A  general  trend  or  indication  of  where  the  fundamental  slope 
Y_-  became  predominant  was  also  of  interest  in  choosing  points  for  the 
least  squares  approximation,  however,  so  all  points  were  initially 
plotted.   These  are  shown  in  Figure  8.   Once  the  apparent  Bragg  cut-off 
break  point  was  fixed  upon,  the  points  making  up  the  two  slopes  were 
submitted  to  separate  least  squares  polynomial  approximations  to  obtain 
a  mean  value  of  slope  and  a  standard  deviation.  By  applying  Chevenet's 
criterion  [9]  for  rejecting  an  experimental  reading,  certain  points 
were  eliminated  for  each  slope.   The  remaining  data  points,  those 
within  the  criterion,  were  submitted  to  a  final  least  squares  fit.   The 
slopes  obtained  in  this  manner  were  entered  in  equation  (3)  to  determine 
the  diffusion  lengths.   The  slopes  y_1   and  diffusion  lengths  L  are 
shown  in  Table  2  for  both  temperature  conditions. 
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TABLE  2 


y,.,  (relaxation  length) 


Temperature 

295°K 
481°K 


Temperature 

295°K 
481°K 


Thermal 


Cold 


0.077  ±  0.013  cm 


-1 


0.026  ±  0.004  cm 


-1 


0.074  ±  0.024  cm 


-1 


0.037  ±  0.008  cm 


-1 


L   (diffusion  length) 


Thermal 


32.72  cm 
45.49  cm 


Cold 

Imaginary 
Imaginary 


B«   FINAL  RESULTS 

As  can  be  seen  from  Table  2  the  diffusion  lengths  are  much  lower 
than  those  normally  found  for  thermal  neutrons  in  graphite.   For 
instance  S.  K.  Trikha,  et  al.  121]  have  calculated  thermal  group 
diffusion  lengths  of  56.75  cm  and  65.17  cm  for  300°K  and  523°K, 
respectively.   Also  they  have  found  a  cold  group  diffusion  length  of 
24.25  cm  and  15.96  cm  for  the  same  temperatures.   In  the  present  case 
the  thermal  group  diffusion  lengths  are  much  lower.   For  the  cold 
group,  L  becomes  imaginary.   That  is,  the  transverse  buckling  term  in 
equation  (3)  is  larger  than  the  relaxation  length  y-i-i  •  DeJuren  and 
Swanson  [5]  also  found  imaginary  diffusion  lengths  in  graphite  piles 
with  sides  less  than  28  inches. 

An  important,  additional  consideration  is  the  fact  that  the  Grade 
HC  graphite  has  many  more  impurities  than  the  AG0T  graphite  normally 
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used  in  experiments  of  this  nature.   Mantell  [17]  lists  the  ash  content 
for  Grade  HC  as  0.30  per  cent  whereas  for  the  AGOT  graphite,  the  ash 
content  is  0.07  per  cent.   Although  no  exact  content  breakdown  listing 
was  found  for  this  particular  grade,  other  manufacturers  of  graphite 
for  similar  applications  list  varying  degrees  of  iron,  vanadium,  sulfur, 
calcium,  titanium,  aluminum  and  boron  as  constituents  of  ash.   The 
presence  of  these  elements  could  have  a  drastic  effect  upon  reducing 
the  diffusion  length  throughout  the  pile. 

Appendix  B  discusses  the  uncertainty  analysis  and  it  is  seen  that 
very  high  inaccuracies  are  present  for  almost  all  data  points,  due  to 
the  low  total  count  received  and  the  relatively  high  background  count. 
It  should  also  be  noted  that  the  steepness  of  the  thermal  group  slope 
makes  it  very  sensitive  to  data  variations. 
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VI.   SUMMARY 

It  has  been  shown  by  theoretical  analysis  that  the  thermal  diffusion 
length  as  illustrated  by  the  relaxation  length  does,  in  fact,  differ 
above  and  below  the  Bragg  cut-off.   Further,  this  difference  decreases 
as  temperature  is  increased,  indicating  that  a  single  straight  line 
(slope),  or  diffusion  length,  is  approached  throughout  the  entire  thermal 
range  at  a  sufficiently  high  temperature  T  ,  where  T  is  greater  than 
1020°K.   The  slopes  listed  in  Table  2  and  shown  in  Figures  8  and  9 
verify  this  analysis.   Although  the  slopes  thus  obtained  for  the  neutron 
group  above  the  Bragg  cut-off  give  lower  than  usual  thermal  group 
diffusion  lengths,  the  trend  is  clearly  brought  out. 

In  the  development  of  equation  (2)  we  have  assumed  exponential  decay 
and  that  the  transverse  leakage  has  no  effect  on  the  total  flux  spectrum. 
However,  as  shown  by  Williams  [22],  for  crystalline  moderators  there  is 
a  maximum  transverse  buckling  which  invalidates  this  assumption. 
Theoretically,  for  the  case  of  graphite,  values  of  buckling  for  sides 
less  than  45  inches  will  lead  to  a  non-exponential  decay  of  the  flux, 
and  its  energy  spectrum  along  the  z  direction  will  change.  No  discrete 
values  of  the  relaxation  length  would  exist  in  this  case.  DeJuren  and 
Swanson  [5]  however  reported  a  diffusion  length  measured  in  a  graphite 
block  with  28  inch  sides.   Only  when  their  sizes  were  reduced  to  20 
inches  did  their  diffusion  lengths  become  imaginary.   As  DeJuren  and 
Swanson  showed  in  their  experiments,  the  transverse  leakage  of  neutrons 
increases  along  the  z  direction  due  to  the  fact  that  cold  neutrons  with 
longer  mean  free  paths  escape  easier  and  therefore  spectral  equilibrium 
is  not  maintained. 
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This  effect  is  apparent  in  the  cold  group  region  of  figures  8  and  9. 
Due  to  the  small  size  of  the  present  pile,  the  spatial  decay  is  no 
longer  exponential  and  equation  (2)  can  no  longer  be  valid. 

The  heating  capacity  of  this  graphite  pile  has  not  been  approached. 
Further  work  in  this  area  at  even  higher  temperatures  should  be  done  to 
verify  this  conclusion,  which  is  based  on  only  one  high  temperature 
test.   It  would  be  recommended  however  that  provisions  for  a  stronger 
neutron  flux  within  the  graphite  be  provided,  in  which  case  the  thermal 
neutron  sensitive  dysprosium  could  still  be  used.   The  possibility 
should  also  be  investigated  of  using  a  water  cooled  BF_  detector  in 
conjunction  with  a  multi-channel  analyzer.   Either  case  will  reduce  the 
uncertainties  presented  here.   Another  important  point  is  that  the 
block  size  must  be  increased  beyond  the  border  line  between  the  imagi- 
nary and  real  diffusion  length  regions.   In  fact,  the  size  -must  be  large 
enough  to  insure  the  existence  of  a  diffusion  length  for  the  cold 
neutron  group. 
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APPENDIX  A  FORMULAE  FOR  DATA  REDUCTION 

In  determining  the  neutron  flux,  use  is  made  of  the  radioactive 
atoms  of  dysprosium  which  are  produced  by  neutron  absorption  within 
the  wire  detector.   The  rate  of  production  R  less  the  rate  of  decay  A, 
called  the  activity,  is  equal  to  the  rate  of  change  of  these  radio- 
active atoms 

dt"  =  R  "  Xn 

where  n  is  the  number  of  radioactive  atoms  and  1      is  the  decay  constant. 
For  steady  flux  cases  R  is  a  constant  and  the  solution  to  the  above 
equation  can  be  written 


A  =  An  =  AgCl  -  e  At)  (7) 


where  A  ,  the  saturation  activity,  is  the  equivalent  of  the  rate  of 
production  R  for  very  long  irradiation  times.   After  removal  from 
irradiation  the  rate  of  decay  will  decrease  exponentially  as 

A  =  An  =  A  e   1 

where  t.  is  the  time  after  removal  and  A.  is  the  activity  at  removal 
time.   If  counting  is  employed  to  determine  the  number  of  radioactive 
atoms  n  existing  over  a  given  period  of  time,  those  counts  C,  less  any 
background  counts  B,  will  be  equal  to  the  difference  between  the  rate 
of  decay  at  the  beginning  t   and  end  t,  of  counting  time,  divided  by 

cl  D 

the  decay  constant  and  multiplied  by  the  efficiency  e  of  the  counting 
equipment. 


41 


n  =  C  -  B  = 


«  ■  ^  ^^a 


x  e 


X(C  -  B)   =  eA   [  e"Ata  -  e"Ub  ] 


(8) 


In  the  case  where  the  detector  wire  is  alternately  removed  and 
inserted  in  a  neutron  flux,  a  recursive  formula  was  derived  for  A.. 
Referring  to  Figure  10,  the  initial  activity  at   t  =  1  will  be  simply 


A,   =  A   (1  -  e  Uil) 
1      s 


(9) 


Between  t  =  1  and   t  =  2  ,  the  value  of  A  will  decrease  an  amount 


A  e   o2  . 


Thus  the  value  of  X  will  equal 


X  =  A1  e"Xto2 


(10) 
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Fig.  10 
Repeated  Exposure  and  Withdrawal  From  Flux 
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Considering  an  exponential  rise  again  at  point  X  and  time  starting  over 
at  t  =  2,  the  value  of  A9  is 

A„   =  X  +  (  A  -  X  )(  1  -  e"Ati2  ) 
2  s 

A0  =  X  e"Ati2  +  A   C  1  -  e~Ati2  )  (11) 

I  s 

Combining  equations  (10)  and  (11) 

a  e~A(to2+ti2)  +  A  (  1  -  e"Ati2  ) 
1  s 

and  from  equation  (9) 


A0     =     A  (   1  -  e"Atil   )    e"A(to2+ti2)     +     A   (  1  -  e"Ui2  ) 
2.  s  s 


A2 


which  can  be  written 


A2     «     AsIl-  e"U"  +   C  1  -  e"Atil  )e"A(to2+ti2)   j 


A  similar  line  of  reasoning  can  be  applied  to  A~ 


A   =  A  I  1  -  e-Xti3  +  A9  e-Uto3+ti3)  j 
3      s  I 


the  recursive  formula  being 

A   -  A  II-  e"Xti3  +  A.  .  e~XCtoJ+tJ45  ]    ,  j=2,3,...   C12) 

3  s  j-1  •  J   • 

Equation  (9)  is  still  valid  for  j  =  1.   From  equations  (8)  and  (12) 
A(C  -  B) 


A 
s 


E   ]        _  e"Ui3  +  A         e"A(t°J+t«)   ]    (   e"Ata  -  e"Atb   )      (13) 


A3 


For  the  case  where  the  flux  is  constant  the  saturated  activity  can 
be  written 


As  -  NT 


J    oa(E)<KE)  dE 


where  N  is  the  total  number  of  target  neutrons  and  a  (E)      is  the 

X  a. 

capture  cross  section  at  energy  E.   Since  the  dysprosium  was  within  two 
per  cent  of  indicating  only  thermal  energy  neutrons,  the  integration 
can  be  carried  out  over  the  thermal  energy  range.   Considering  a 
Maxwellian  flux  and  assuming  a  1/v  behavior  for  dysprosium,  the 
saturation  activity  can  be  written 


As  "  Va«  Sr  >  1/2  •   (1A> 


where  T  is  293.16°K,  T  is  the  medium  temperature,  and  a  ■  920  b. 
is  the  capture  cross  section  recommended  in  Ref.  7.   It  is  based  on  the 
weighted  cross  sections  for  Dy  isotopes  by  their  abundance  in  natural 
dysprosium.   Equations  (13)  and  (14)  now  give,  for  the  thermal  neutron 
flux 

♦  -  _ii£^ (15) 

ttT 
eA.    C  e"Ata  -  e"Atb   )   0&   (  ^  )    1/2  NT 

N  ,  the  total  number  of  target  neutrons,  can  be  found  by 

NT  =  V  (  m/V   '   A°  )  v 
1  A 

where  V  is  the  wire  volume,  WT  is  the  wire  weight,  A  is  Avogadro's 
number,  A  is  the  average  atomic  mass,  and  v  is  the  fraction  by  weight 
of  dysprosium  in  the  Dy-Al  wire. 
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APPENDIX  B   UNCERTAINTY  ANALYSIS 

As  stated  earlier  the  low  thermal  flux  resulted  in  a  relatively  low 
count  rate  compared  to  the  background  count.   An  analysis  was  conducted 
on  each  of  the  data  points  to  determine  the  accuracy  of  the  resulting 
flux. 

Considering  a  total  of  C  counts  as  the  mean  value,  its  standard 

1/2 

deviation  can  be  written  as  C    ,  and  the  total  count  would  be  reported 

1/2 
as   C  ±  C    .   The  same  criterion  is  applied  to  the  background  count  B, 

In  obtaining  the  difference,  C  -  B,  over  the  same  time  interval,  the 

value  with  its  standard  deviation  is 


(C  -  B)   ±   (ac2  +  aB2)  1/2 


where  the  standard  deviations 

are 

°C  = 

c1/2 

and 

aB  = 

B1/2 

This  value  is  converted  into  a  percentage  error  by 

t      2  j    2,  I/2 

(C  -  B)  ±  — — 100% 

(C  -  B) 

Since  this  was  the  major  source  of  uncertainty  in  determining  the  flux, 
the  same  percentage  error  applies  to  the  thermal  flux  obtained  for  each 
data  point.   These  uncertainties  in  terms  of  percentage  are  tabulated 
in  Tables  3  and  4  at  the  end  of  the  Appendix. 

When  the  data  points  were  plotted  there  were  some  obviously  erratic 
points.   It  seemed  appropriate  to  disregard  some  of  the  more  erroneous 
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ones.   Therefore  a  least  squares  approximation,  based  on  equation  (6), 

was  obtained  on  the  computer.   The  standard  deviation  of  the  total  of 

all  data  points  was  also  obtained  from  the  best  fit  straight  line. 

Using  Chevenet's  criterion  for  rejection  of  outlying  data,  a  ratio 

d. 

( — )  of  data  deviation  from  its  mean  value  to  the  standard  deviation 
o 

was  obtained  for  each  position.   Those  points  giving  ratios  outside  the 
criterion  were  rejected  and  those  less  than  the  criterion  were  resubmit- 
ted for  a  second  least  squares  approximation.   It  is  this  second 
approximation  for  slope  that  is  used  for  the  final  results  and  discussion. 
The  ratios  obtained  as  above  are  also  listed  in  Tables  3  and  4. 

The  high  temperature  uncertainty  was  determined  as  follows,  utilizing 
the  voltage  readings  in  all  calculations.   When  the  final  value  was 
obtained,  the  conversion  to  degrees  centigrade  was  applied.   At  each 
of  the  levels,  upper  and  lower, 

1=1 

T    1       S     i~  n2  T   1/2 

1=1 

where  v .  are  each  of  the  five  thermocouple  readings  at  one  recording 
time  and  N  is  equal  to  the  number  of  thermocouples.  Then  the  mean  and 
deviation  for  the  pile  at  each  recording  instance  is  found 

v  =  \   (uu  +  vt) 

1  ,  2     2.  1/2 
av  "  2  (au  +  °l  ) 

where  the  subscripts  u  and  I   refer  to  upper  and  lower  levels, 

46 


respectively.   The  overall  temperature  value  and  its  standard  deviation 
is  then  obtained  from  the  mean  values  of  each  recorded  run: 


i  ;  - 

Vm     =     N       *       Vi 

1=1 


a  E       a     2 

m  .  ,       v. 

i=l         l 


Using  the  above  procedure  a  high  temperature  of   207.88  ±  A.56°C  was 
obtained. 
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TABLE  3 


Temperature:      295°K 


Height  above 
base   (cm) 

Flux 
(n/cm^sec) 

Error 
(per  cent) 

a 

Chevenet ' s 
Criterion 

35.21 

2354.8 

8.60 

0.05 

1.73 

37.65 

1926.5 

11.61 

0.57 

1.73 

40.11 

1686.2 

11.25 

0.14 

1.73 

42.57 

1311.7 

17.30 

1.25 

1.73 

45.02 

1523.3 

12.62 

4.51 

1.73 

47.53 

849.6 

25.35 

2.79 

1.73 

49.99 

896.9 

20.25 

2.93 

2.13 

52.43 

978.9 

21.30 

0.14 

2.13 

54.85 

761.4 

22.75 

4.12 

2.13 

57.29 

798.4 

24.95 

1.23 

2.13 

59.75 

692.7 

21.85 

1.64 

2.13 

62.21 

752.9 

23.50 

0.69 

2.13 

64.63 

557.5 

27.20 

3.19 

2.13 

67.03 

523.9 

32.90 

3.52 

2.13 

69.47 

492.4 

30.90 

3.69 

2.13 

71.91 

748.2 

24.30 

3.82 

2.13 

74.35 

436.3 

34.70 

9.26 

2.13 

76.81 

563.1 

31.50 

4.01 

2.13 

79.25 

439.8 

35.20 

4.59 

2.13 

81.67 

389.9 

45.90 

0.64 

2.13 
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TABLE  4 


Temperature:   481  °K 


Height  above      Flux        Error 
base   (cm)     (n/cm^sec)     (per  cent) 


32.75  5518.9  10.35 

35.21  8612.5  13.55 

37.65  2204.9  20.23 

50.11  1713.0  12.17 

42.57  2169.6  12.17 

45.02  2502.8  15.74 
47.53  1671.0  19.00 
49.99  1547.5  24.40 
52.43  1204.2  30.50 
54.85  3291.3  31.55 
57.29  890.6  30.90 
59.75  1055.6  26.50 
62.21  1050.0  39.75 
64.63  802.9  36.30 

67.03  740.4  30.45 
69.47  1022.5  32.10 
71.91  828.8  152.50 
74.35  646.8  115.00 
76.81  222.1  149.50 
79.25  466.9  54.90 
81.67  (neg)  74.70 


X 

a 

Chevenet  *  s 
Criterion 

0.32 

1.90 

4.93 

1.90 

3.48 

1.90 

3.87 

1.90 

0.78 

1.90 

1.66 

1.90 

0.19 

1.90 

1.04 

1.90 

1.76 

2.10 

3.01 

2.10 

7.79 

2.10 

3.59 

2.10 

0.79 

2.10 

0.32 

2.10 

1.11 

2.10 

0.74 

2.10 

3.52 

2.10 

2.59 

2.10 

1.43 

2.10 

6.50 

2.10 

1.79 

2.10 
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